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We apply chiral quark model with momentum dependent quark mass 
to two kinds of nonperturbative objects. These are: photon Distribution 
Amplitudes which we calculate up to twist-4 in tensor, vector and ax- 
ial channels and pion-photon Transition Distribution Amplitudes together 
with related form factors. Where possible we compare our results with 
experimental data. 

1. Introduction 

One of the biggest problems in particle physics is description of hadrons 
in terms of the fundamental degrees of freedom — quarks and gluons. This 
is in fact a non-perturbative problem and usually is formulated in terms of 
various distribution functions, which appear in QCD factorization theorems. 
The most famous example are Parton Distribution Functions which can be 
measured in the inclusive lepton-hadron deep inelastic processes. They are, 
however, one dimensional distributions only. Therefore although they are 
by now sufficient for description of various processes at high energies, they 
simply give only limited information on the structure of hadrons. 

On the other hand one can study also hard exclusive processes, such 
as deeply virtual Compton scattering for instance. Then the factorization 
theorem states, in great simplicity, that the amplitude is given by the con- 
volution 

M = (soft) ® (hard) , (1) 

where hard is the part that can be calculated in perturbative QCD, while the 
soft part is of non-perturbative nature. In the following we will be mainly 
concentrated on the soft part. Although difficult to access experimentally, 
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they can be obtained either by lattice calculations or — as we shall see 
- they can be estimated from theoretically justified low energy effective 
models. 

The soft part parametrizes hadronic matrix elements of certain non- 
local operators on the light-cone. The simplest objects of this kind are Dis- 
tribution Amplitudes (DA) which correspond to hadron-to-vacuum matrix 
elements of bi-local quark operators on the light-cone. In case of the leading 
twist DAs, i.e. the ones giving main contribution to the amplitude, they 
describe (in the infinite momentum frame) the probability for a composite 
particle to dissociate into its constituents with given longitudinal momentum 
fractions. Distribution Amplitudes have been successfully used in theoretical 
description of hadronic form factors [U [2] for many years. However recent 
BaBar data for pion-photon transition form factor shows, that probably the 
standard factorization formulae do not apply O 0] . We shall come back to 
the BaBar data in Section [H 

More general class of soft objects are Generalized Parton Distribu- 
tions (GPD). They correspond to non-diagonal in momenta matrix elements, 
therefore they describe also the distributions of transverse momenta of the 
partons inside the hadron. GPDs appear in description of deeply virtual 
Compton scattering for instance, which is recently the subject of intensive 
theoretical and experimental studies. For a review of this issue see e.g. [7]. 

One can still define more general class of the objects than GPDs — so 
called Transition Distribution Amplitudes (TDA). They parametrize matrix 
elements which are non-diagonal in momenta and in physical states. Such a 
family of objects was introduced for the first time in [13J. We shall discuss 
this class further in Section |4) 

As already remarked above there is very little experimental data concern- 
ing the soft part of ([1]). On the other hand it can be studied in effective 
models. This is however nontrivial not only because of complex non-local 
interactions at low energies. Even bigger problem is that in general effective 
models do not inherit all symmetries of the underlying theory. Soft objects 
considered here appear in the framework of QCD, therefore they should 
posses several important properties, for example Lorentz and gauge invari- 
ance. They should also correctly reproduce quantum anomalies. The task 
to cope with all the constrains in the effective models is therefore nontrivial. 
We shall come back to this point in Section [2] 

2. Non-local chiral quark model 

Let us consider the scattering process involving the simplest possible 
hadronic state — the pion. On one hand it is a bound state of quark- 
anti-quark pair and the Goldstone boson of spontaneously broken chiral 
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symmetry on the other. Before we proceed let us briefly recall these very 
important aspects of QCD. 

Here and in the following we assume only two quarks u and d which 
are massless, i.e. m u = = 0. Then the Lagrangian of QCD is in- 
variant under separate rotations of left- and right-handed spinors, that is 
the symmetry group is SU (2) R ® SU (2) L (the chiral symmetry group). It 
is generated by the chiral charges satisfying SU (2) commutation relations 

Qh R = / ( ^' x R TSTtV^R ( x )i where r a are Pauli matrices, tp denote 
iso-doublets. One can also define the combination of chiral fields transform- 
ing as vector and axial-vector. Then the corresponding combination of L, R 
charges Q a = Q R + Q£ and Q% = Q R - Q£ generate the SU (2) v ® SU (2) A 
group3- The most direct consequence of this symmetry would be a de- 
generacy of the states with different parity. However, such a behavior is 
not seen in the hadronic spectrum — on the contrary, we observe huge 
mass differences between parity partners. The most natural way of solving 
this discrepancy, is to postulate that although theory is chirally invariant, 
the vacuum state is not. This phenomenon is known as spontaneous chi- 
ral symmetry breaking (SxSB). According to Goldstone theorem we should 
then observe a triplet of massless pseudo-scalar particles — the Goldstone 
bosons. Indeed they can be apparently identified with pions (7r + , 7r°, 7T~), 
which are very light (m n ~ 140 MeV) in comparison to other hadrons (e.g. 
^proton ~ 1 GeV). Non-zero pion mass can be explained by finite (although 
small) current masses of u, d quarks, which explicitly break chiral symmetry 
from the very beginning. 

Another important aspect of SxSB is the existence of the quark conden- 
sates, i.e. the quantities 

(0 \qq\ 0) = (qq) = (q R q L ) + (q L q R ) , (2) 

where q denotes either u or d quark field. It can be easily seen that the 
nonzero value of the quark condensate breaks chiral symmetry of the vac- 
uum. Consider the commutator 

Q^7 5 tV] =-5 ab ^ (3) 

and its vacuum expectation value. If the right hand side is nonzero it implies 
that 

Of 10)^0, (4) 

what is exactly the SxSB condition. Therefore the quark condensate can be 
viewed as an order parameter measuring the breakdown of chiral symme- 
try. Phenomenological value of the quark condensate is quite large (qq) ~ 



There is also similar global symmetry acting on the whole doublet. The axial sym- 
metry U (1) A is however broken due to quantum anomaly. 
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(-250 MeV) 3 (at renormalization scale about lGeV). Let us notice next 
that in QCD (qq) is represented by a closed quark loop, i.e. it is propor- 
tional to the trace of fermionic propagator (qq) ~ TrS(x,x), where the 
trace is over Dirac and color indices. However, if this quantity is non-zero 
there must be a non-slash term in the propagator — the mass term. This 
dynamically (due to SxSB) generated mass is often referred to as constituent 
quark mass. Notice that the quark condensate is a purely non-perturbative 
quantity, since it is impossible to generate non-slash quark self energy by 
interactions of vector bosons. Rather it must be created by some kind of 
a scalar interactions. We shall come back to the issue of quark condensate 
later in Section [3l 

Let us now switch to description of the interactions between the pions 
and quarks at low energies. It is clear from the above that such a model must 
incorporate SxSB. It is convenient to start discussion by recalling the famous 
Nambu — Jona-Lasinio (NJL) model. It is an effective theory of quarks with 
four fermion couplings, appearing due to integrating out the gluonic degrees 
of freedom from the QCD action. In the standard NJL model couplings with 
more fermions are neglected. The Lagrange density for the simplest version 
of the NJL model reads 



(5) 



where r = (r 1 ,r 2 ,r 3 ) are Pauli matrices and G is coupling constant. It can 
be checked using some algebra and the relation e - *^" 7 ") 75 = (cos \a\ — 275a • r sin 
where a 1 = \a\a l that Lagrangian ([5]) is indeed chirally invariant. The most 
important feature of NJL model is that it incorporates the mechanism lead- 
ing to SxSB. One way to see this is to solve the corresponding lowest order 
Dyson-Schwinger equation for the quark propagator. Denoting quark self- 
energy by £ (p) = M one obtains the following consistency condition (so 
called gap equation) 



It has two solutions: M = (for massless quarks) and The latter 

corresponds to the constituent quark mass which generates non-zero quark 
condensate breaking the chiral symmetry of the vacuum. Notice that the 
integral in ([6]) requires regularization. We shall discuss this later in this 
section. For a review of NJL model see [6] for example. 

Mesons can be easily introduced into the just described theory as aux- 
iliary fields a and 7r a — this can be done formally in the path integral 
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formalism and is called bosonization procedure. The new Lagrange density 
reads 




where g 2 = fi 2 G. Notice that the fields a, ir a are truly auxiliary - there 
are no corresponding kinetic terms, moreover they are composed fields what 
can be immediately seen using equations of motion. NJL Lagrangian in the 
form ([7]) can also be used to show that the ground state which minimizes 
the energy is populated by the scalar quark condensate. 

One can also look at the appearance of mesonic fields from a slightly 
different point of view. Consider the following effective Lagrange density 

C = $(i$- M)ip, (8) 

which leads to Dirac equation for the quark with constituent quark mass M. 
However ([8]) is obviously not chirally invariant. In order to fix this deficiency 
one has to introduce additional fields in the form 

LH5 ( x ) = e ^-n(x) l5 _ 1 + ( s ) + _ (g) 

where F n ~ 93MeV is the pion weak decay constant, and couple them to 
quarks, 

C = ${i$- MW 5 )il). (10) 

Then the axial transformations of quark fields can be absorbed by pion fields 
ir a . The Lagrange density (|10p is a starting point for our further consider- 
ations and represents the simplest local chiral quark model. It describes 
quarks having dynamically generated constituent mass M and interacting 
with the external pion fields. 

The effective theories just described are a non-renormalizable ones. The 
regularization introduced in order to make the loop integrals finite cannot 
be removed at the very end of the calculations and the observables depend 
on its actual form. Moreover it is somehow (but not straightforwardly) 
related to the domain of applicability of the model. There are many ways 
of regularizing the loop integrals. One could use for example simple four- 
momentum cutoff or Pauli-Villars regularization. However, the point is that 
the regularization scheme should respect all symmetries of the underlying 
theory, i.e. QCD. This is extremely important especially in the case of soft 
matrix elements as stated in Section [TJ Therefore four-momentum cutoff is 
excluded in the first place since it violates Lorentz invariance. Also very 
often used Pauli-Villars regularization is not the best method, because in 
order to get the results consistent with QCD one must keep it finite in some 
diagrams and remove in others (connected with anomalous processes). On 
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the other hand notice that in reality the constituent quark mass should not 
be a constant — it should vanish for large quark momenta due to asymptotic 
freedom. Therefore in the following we assume that 

M = M (k) = M F 2 (k) , where F (ft) — ► 0, F (0) = 1. (11) 

The constituent quark mass at zero momenta M$ is chosen to be about 
M ~ 350 MeV. 

The interaction part of the effective action corresponding to (|10p with 
assumption can be written in momentum space as 



Sint = M J ^^L^k)F(k)U^(k-l)F(l)m- (12) 

The explicit shape of F (ft) cannot be obtained from the gap equation itself. 
However the action (|12|) was actually obtained in the instanton model of the 
QCD vacuum, together with the expression for F (ft) = Fi nst (ft). Unfortu- 
nately -Finst (ft) turns out to be a highly non-trivial function of Euclidean 
momenta [8]. Therefore instead of -Fi ns t (ft) we shall use the following simple 
formula in Minkowski space 



F<t) = l^f?17) • (13 > 



which reproduces -Fi ns t quite well when continued to Euclidean space. The 
parameter n is responsible for the actual shape of F(k), therefore we can 
investigate the sensitivity of calculated quantities to the form of the cutoff 
function. The cutoff parameter A n is adjusted in such a way that pion decay 
constant F n given by the formula [21 J 



^ _ N c f°° , l2l 2 M 2 fas) - k l M ( fc E) M' (k E ) + k%M' (k E ) 



n Wo (fel + M2(fe E )) 2 

is equal to the experimental value. In the above equation k E corresponds 
to Euclidean momentum, while the prime to differentiation with respect to 
k E . For example for constituent quark mass Mq = 350 MeV and n = 1 (|14[) 
gives Ai = 836 MeV. 

Although momentum dependent quark mass seems to be the most nat- 
ural regulator it introduces a serious difficulty Namely Ward-Takahashi 
identities are not satisfied in such a model. It can be most easily seen by 
considering a divergence (in momentum space) of the vector current and 
applying Dirac equation. This violation turns out to be not very large, it 
can however spoil some important properties of soft matrix elements, like 
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correct normalization for example. In order to fix this problem the standard 
vector 7^ and axial 7^75 vertices have to be modified by adding new non- 
local terms. The problem is, however, that such a modification is not unique 
|21| (Ward identities fix only longitudinal part of the vertices) . In this work 
we use the following modified vector and axial vertices: 

(k,p) = - -^-liL (M {k) _ M (p)) , (15) 



It (M = 7^75 + 7 TT2 ( M ( k ) + M (?)) 75, (16) 



(p-kf 

reproducing Word-Takahashi identities. The vector vertex does not intro- 
duce additional singularities, while the axial one has a pole corresponding to 
the massless pion as it should be [20] . Let us remark that eq. (|14p expressing 
F n is determined unambigously since it involves only the derivative of the 
axial current. 

At the end of this section we remark that any non-local (i.e. with mo- 
mentum dependent constituent quark mass M (k)) chiral quark model is 
determined by specifying both the M (k) and the precise form of all ver- 
tices. 



3. Photon Distribution Amplitudes 

As already remarked in the Introduction the simplest soft objects are 
Distribution Amplitudes. In this section we present how the non-local chiral 
quark model can be applied to this class. However, instead of considering the 
hadronic DA we shall discuss less known photon DA. This is possible due to 
the fact that — besides standard perturbative part — photons possess also 
hadronic component. This fact is very well known from photoproduction 
processes, where photon structure function has to be taken into account (so 
called resolved photoproduction). 

Photon DAs appear for example in the description of vector mesons 
radiative decays. To be more specific consider for instance the process 
D°*(q + p) ->■ D° (p) +j(q) (Fig. [3k.). Using OPE one can then write 
the amplitude as products of factors that are divergent on the light-cone 
and finite photon-to- vacuum matrix elements. The latter can be identified 
with the photon DA as we shall see. 

Before we give the more precise definition of the photon DA we should re- 
call usefull kinematical variables. One defines two null vectors n = (1, 0, 0, — 1) 
and h = (1, 0, 0, 1). Then any four- vector can be decomposed into "plus", 
"minus" and transverse components 

= v+— + v-— + v£. (17) 
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a) b) 




Fig. 1. a) Bag diagram for the radiative D°* vector meson decay. The lower blob 
corresponds to photon Distribution Amplitude, b) Simple quark loop correspond- 
ing to photon DA in the quark model. Double external line represents bilocal quark 
operator on the light-cone. 

Photon Distribution Amplitude is denned as a Fourier transform of the 
photon-to-vacuum matrix element of the non-local quark operator on the 
light-cone. In general this can be written as 

J ^ e ^ u -^ p+ <0 p(An) (-An)| 7 (P)> ~ F a (P 2 ) 

X {Otwist-2 <Po 1St - 2 («, P 2 ) + Otwist-3 <t% iSt " 3 («> + • • •} ( 18 ) 

where O = {c^ u , 7^, 7^75} corresponds to different tensor nature of bilocal 
operators, O tw ist-2, Otwist-3, • • • denote apropriate tensor structures which 
are multiplied by photon DA 4>o 01 given kinematical twist. Notice that we 
do not assume that the photon is on-shell. Then the decay constants Fq 
depend on photon virtuality P 2 and become a kind of "form factors" — we 
shall use this terminology in the following. For more precise definitions of 
the photon DAs refer to P HP) [17] . 

Using the non-local chiral quark model described in Section [2] we calcu- 
lated photon DAs up to twist-4 in tensor, vector and axial channels |12| in 
one loop approximation (Fig. [3b). Our results are analytical up to the solu- 
tion of a certain polynomial equation. Let us briefly summarize our results. 
For real photon the leading DA is the twist-2 tensor amplitude <f) a v = (J)t 
(Fig. [3^) corresponding to a^ u structure. We find it is almost flat and non- 
vanishing in the end-points. Also the sensitivity to the n parameter, i.e. to 
the shape of F (k) is rather small. In the vector channel one has to subtract 
the infinite, perturbative part when calculating the corresponding matrix 
element. Then we find in particular that leading twist vector DA vanishes 
in the end-points. It can be easily shown on general grounds that the vector 
"form factor" P 7 m (P 2 ) = Fy (P 2 ) should be zero for the real photon. This 
property is maitained in our model only when we use modified vector vertex 
(Fig. [3b), as described in Section [2l Higher twist amplitudes turn out to 
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Fig. 2. a) Twist-2 tensor photon Distribution Amplitude for n = 1 and several 
values of the constituent quark mass Mq. b) Vector form factor for Mq = 350 MeV, 
n = 1, calculated using naive vector vertex 7 M (dashed line) and the modified one 
T' 1 (solid). Notice that the modified vertex assures that Fy vanishes for real photon 
as required by QED. 



be rather strongly model dependend. Moreover some of them contain Dirac 
delta functions in the end points, they should be therefore viewed rather as 
the generalized functions. Similar calculation was previously done in Ref. 
|10) and differs from ours in some points. 

The left hand side of the definition (|18p is dimensionfull, therefore we 
should have several quantities that set up the characteristic mass scale for 
photon DAs. Among others, it is a quark condensate, already disscussed in 
Section [2j The non-local chiral quark model with f)13[) allows to obtain the 
following "analytical" expression for the quark condensate 

at M 2 A 2 4 «+! 

m = - \^ n E f^ n o- + ln c 1 + > ( 19 ) 

2=1 

where the complex numbers r]i are numerical solutions to the equation 
z 4n+i +z 4n _ ( Mo/An ) 2 = 0, while are defined as /, = Ut^f ~ 

For example for M = 350 MeV and n = 1 we get (qq) = (-253 MeV) 3 . It 
turns out that in general the values of (qq) rather strongly depend on model 
parameters. 

4. Pion-photon Transition Distribution Amplitudes 

In this section we switch to more involved applications of the non-local 
chiral quark model. Transition Distribution Amplitudes, apart from be- 
ing interesting on their own, can serve as a demanding testing ground for 
the model. The reason is that they involve diagrams responsible for axial 
anomaly. We shall come back to this point later in this section. 
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Fig. 3. a) The bag diagram for the process 7r + 7r~ — > 7*7. The lower bag represents 
Transition Distribution Amplitude while the upper corresponds to the hard process, 
b) The quark loop corresponding to TDA, the bilocal operator is assumed to "live" 
on the light-cone. In order to recover correct normalization both vertices have to 
be non-local. 

Transition Distribution Amplitudes were originally introduced in order 
to describe hadron-antihadron annihilation into two photons, i.e. the pro- 
cess HH — > 7*7 or backward virtual Compton scattering j*H — > |13j . 
The amplitudes for these processes can be described in QCD analogously 
to the reactions HH — > 7* and j*H —> H respectively with the restriction 
that Distribution Amplitudes for H should be replaced by a new object - 
Transition Distribution Amplitudes (Fig. (3^i). First estimates were done in 

Refs. [HUH USE]. 

Before we give the general definition of TDAs we should define relevant 
kinematics. We consider pion with momentum P^ transforming into the 
photon with momentum Pi^ . We define the momentum transfer as q^ = 
P% — Pi and the momentum transfer squared t = q 2 which is assumed to 
be small. Using the average momentum p^ = \ {P\ + P%) we define so 
called skewedness £ = —q + /2p + , which is a standard variable in the GPDs 
formalism. We consider chiral limit and real photons, i.e. P 2 = P 2 = 0. 

The general definition of leading twist TDAs can be written as 

/ <T ( p 2,£) I? (An) (-An) | vr+ (P 1 )) = O twist _ 2 D (X, £, t)+. 

where in practice O = {7^,7^75}- Dots stand for the other terms that 
can appear and are not related to TDA under consideration. For example 
in the axial channel, i.e. for O = 7^75, pion DA accompanied by massles 
pole appears on the right hand side. This reflects the fact that the axial 
current couples to a pion directly. In the following we denote vector TDA as 
V (X, £, t) (i.e. for O = 7^) and the axial TDA as A (X, f , t) (for O = 7^75). 
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a) 



b) 




Fig. 4. a) Vector Transition Distribution Amplitude for M = 350 MeV, n = 1, 
t = —0.1 GeV 2 and £ = 0.5. Solid line corresponds to the full non-local model with 
non-local vertices and is a sum of the dashed line and the dotted. Dash-dotted 
line was obtained in local model, i.e. with M (k) = M. b) The same for the 
axial TDA. Here the addition coming from the non-local part of the vertices gives 
negative contribution. 

There is very important property that TDAs should posses, namely so 
called polynomiality 



which follows simply from Lorentz invariance. In principle the zeroth mo- 
ment is related to the corresponding form factor. Second very important 
constraint is the normalization of the vector TDA, which is fixed by the 
axial anomaly 



Above condition is model independent and can be derived using Ward- 
Takahashi identities that relate the two-photon matrix elements of the axial 
and pseudoscalar currents. The latter can be then identified with our matrix 
element (|20p with O = 7^ . There is no similar normalization condition for 
axial TDA. However, in the local models, i.e. with M (k) = M it turns out 
that 



In the quark model, calculation of the TDAs reduces to performing quark 
loop shown in Fig. [3jo. We present a typical results in Fig. |4]|18j. Notice 
first that curves obtained in non-local model are much more smooth then 




(21) 






(23) 
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the ones obtained in the local model. Next, we find that the normalization 
condition (122 p is recovered only when light-cone bilocal current in ([20]) is 
also modified according to (|15p . At the same time the normalization of 
the axial TDA is much lower than (|23j) . This result is important because 
zeroth moments of vector and axial TDAs are directly related to the vector 
and axial form factors which can be estimated experimentally. To be more 
precise the relation is 



dX 



v(x,c,t) 

A(X,£,t) 



2V2F W 



(t) 



(24) 



-i 



where the superscript \ denotes that these quantities are defined in the 
chiral limit. The experimental values for t = are (PDG) 



p exp 



A 



(0) 
(0) = 



(25) 
(26) 

(F A (0)/F y (0)) exp = 0.7+^. (27) 
On the other hand the normalization ([22[) gives (model independent) 

Fy (0) « 0.027, (28) 



= 0.017 ±0.008, 

0.0115 ±0.0005, 
0.7 ' " 



what overshoots (|25p more than one standard deviation. The results for 
axial form factor are model dependent. For reasonable model parameters 
we obtain: 



M [MeV] 


n 


F*(0) 


n (o) m (o) 


225 


1 


0.0217 


0.80 


350 


1 


0.0168 


0.62 


350 


5 


0.0163 


0.60 


400 


1 


0.0161 


0.60 


400 


5 


0.0152 


0.56 



We see that indeed the assumption of non-locality ([lip lowers the value of 
factor i ? 7r7 via the relation 



towards the experimental data. 

Moreover Fy is directly related to so called pion-photon transition form 



F^ (t) = V2F* (t) . 



(29) 



This quantity describes the pion decay tt° — > 7*7 process and was measured 
by CLEO [23], CELLO [22] and recently by BaBar [2l] collaborations. We 
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Fig. 5. The experimental data for pion-photon transition form factor i ? 7r7 times 
momentum transfer t. The shaded area represents the predictions from the non- 
local chiral quark model predictions for sensible model parameters. We get the 
best description of the low momentum tranfer data for Mo ~ 300 MeV. 

compare our predictions to the experimental ones in Fig. [5j There is however 
important remark in order. Notice, that by definition TDAs are sensible only 
for small momentum transfers t, the precise range of application is however 
not known. Therefore, as an example we have chosen arbitrarily the range 
of — 8GeV 2 . It is worth noting at this point that the new BaBar data 
are in disagreement with the standard QCD factorization formula, as it was 
already remarked in Introduction. In QCD, the pion-photon form factor can 
be described using pion DA and some perturbatively calculable factor, which 
leads to the certain asymptotic form. New BaBar data cover the range of 
— 40 GeV 2 (in Fig. [5] we retained only the relevant low momentum data) 
and cross the asymptotic line already at about 10 GeV 2 . One way to resolve 
this discrepancy is to note that pion DA which vanishes at the end-points 
was assumed in the standard factorization formula. In Refs. 013] the autors 
study the flat pion DA in order to describe the new BaBar data (but see 
also [5])- 



Let us briefly summarize our presentation. In the begining we recalled 
chiral quark models, starting from widely known Nambu — Jona-Lasinio model 
We argued that spontaneous chiral symmetry breaking is their main ingre- 
dient. We also showed that they lead to nonzero quark condensates and in 
turn to dynamically generated constituent quark mass, which in general can 
depend on momentum. Next we used a simple ansatz for this dependence 
and applied the model to two low-energy objects: photon Distribution Am- 
plitude and photon-pion Transition Distribution Amplitude. We find that 



5. Summary 
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they fulfil most symmetries required by QCD, provided we modify the vector 
and axial vertices in such a way that relevant currents are conserved. We 
find also that form factors which are calculated using Transition Distribution 
Amplitudes are realistic when compared to the experimental data. 

At the end we draw attention to important issues which was not cov- 
ered by this presentation. First of all, in QCD all the low-energy quantities 
depend on some factorization scale \x and are a subject for corresponding 
QCD evolution. On the other hand, within effective models they are ob- 
tained at some fixed fx, which is in fact unknown (although can be roughly 
estimated) . Therefore before one makes a real use of them the evolution has 
to be applied. The second remark is rather a technical one and concerns the 
cutoff A n parameter in (|13p . One should not confuse it with the scale /i of 
the model, as discussed in |19] , 
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